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Abstract. 

Non-time reversible phonobreathers are non-linear waves that can transport energy in coupled oscillator 
chains by means of a phase-torsion mechanism. In this paper, the stability properties of these structures 
have been considered. It has been performed an analytical study for low-coupling solutions based upon the 
so called multibreather stability theorem previously developed by some of the authors [Physica D 180 235]. 
A numerical analysis confirms the analytical predictions and gives a detailed picture of the existence and 
stability properties for arbitrary frequency and coupling. 
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1. Introduction 

One of the subjects where a great deal of attention has been focused in the last two decades is the dynamics 
of nonlinear lattices. Intrinsic localized modes or discrete breathers is one of the most outstanding structures 
that arise in those lattices [1, 2]. Discrete breathers are periodic and localized solutions whose existence is 
allowed by the interplay between discreteness and nonlinearity. The existence of those structures in networks 
of anharmonic oscillators (also known as Klein-Gordon lattices) was firstly proven at 1994 by MacKay and 
Aubry [3]. They established that discrete breathers can be continued from the anti-continuous limit (i.e. the 
limit where there is no coupling between the oscillators) to finite coupling as long as no integer multiple of 
the breather frequency resonates with the linear modes in the phonon band. 

This theorem also demonstrates the existence of breathers with more than one excited sites, dubbed as 
multibreathers. These multibreathers receive special names in some cases; for instance, breathers with all 
their sites excited receive the name of phonobreathers. If all but one of the sites are excited the multibreather 
is called a dark breather [4], in analogy to the dark solitons existing in the Discrete Nonlinear Schrodinger 
(DNLS) equation [5, 6]. 

In order to observe experimentally discrete breathers, they have to be stable. The stability of one-site 
breathers (i.e. breathers with only one excited site) was firstly proven in [7, 8]. It was not until 2003 where 
some of the authors of the present paper developed a method for determining the stability or instability of 
many kinds of multibreathers [9, 10] based on the Aubry's band theory developed in Ref. [7]. An alternative 
approach, introduced in [11], is based in previous work by MacKay et al [12, 13, 14]. Recently, the equivalence 
between both approaches has been proven [15]. 

Most of theoretical work related to discrete breathers has been carried on time-reversible solutions of 
Hamiltonian lattices. Non-time-reversible solutions have mostly been considered in dissipative lattices [16, 
17, 18, 19], which have many experimental applications in chains of coupled pendula [20], micromechanical 
arrays [21], transmission lines [22, 23] and Josephson junction arrays [24, 25]. 

However, there have been very few approaches to non-time-reversible breathers or multibreathers in 
Hamiltonian lattices. The existence of one of such solutions was firstly proven in the pioneer work of 
MacKay and Aubry [3]. It was restricted to phonobreathers where there is a constant phase difference 
between adjacent sites. These solutions can transport energy by means of phase torsion and can be viewed 
as nonlinear phonons or phasons. Later on, Aubry [7] demonstrated the existence of a generic non-time- 
reversible breather, independent on the excited sites. The proof was also extended to inhomogeneous lattices 
(e.g. with vacancies) and vortices. It was also proven that phonobreathers with phase torsion generate a 
stationary flux, and that time-reversible phonobreathers cannot transport energy. Finally, Cretegny and 
Aubry numerically demonstrate [26, 27] the existence of non-time-reversible phonobreathers in homogenous 
and inhomogeneous (with vacancies) 1-D lattices and vortices and breather "rivers" (i.e. percolating clusters 
of breathers connecting two boundaries submitted to phase torsion) in 2D square lattices. They also sketch 
some properties of the flux and its stability in ID lattices. Further applications of this theory can be found 
in many publications on discrete vortices in DNLS [28, 29, 30, 31, 32] and Klein-Gordon [15, 33] lattices. 

The aim of this paper is twofold. On the one hand, we will prove by making use of the techniques 
developed in [9, 10], a result mentioned in passing on Ref. [26]: the stability of phonobreathers depends on 
the phase difference between adjacent sites and the phonobreather frequency. On the other hand, we will 
make an analysis of the stability and flux dependence on the system parameters for two kind of potentials 
in Klein-Gordon lattices, checking the validity of the analytical predictions. 

The paper is organized as follows: in Section 2 we introduce the model equations; section 3 deals with 
the analytical results regarding the stability for low coupling; in order to check these results, we make an 
exhaustive numerical analysis of the existence and stability of phonobreathers at finite coupling in Section 
4; finally, we present our conclusions and some possible extension of this work in Section 5. 
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2. Model setup 

2.1. Dynamical equations and energy flux 

We a consider a Klein-Gordon chain of oscillators with nearest-neighbours harmonic coupling. The dynamical 
equations are of the form: 

u n + V'(u n ) + e(2u n - u n+ i - u-i) = n=l,...,N (1) 

where the variables u n are the displacements with respect to the equilibrium positions, V(u n ) is the on-site 
potential, TV is the number of oscillators, and e > is the coupling constant. 

We consider two paradigmatic cases of on-site potentials: Morse (soft) potential 

V(u) = \[ C M-u)-l] 2 (2) 

and (j) 4 (quartic) hard potential 

V(u) = \u* + \u A (3) 

We look for non-time- reversible solutions of the dynamical equations (1) with all the particles excited, 
so that the difference between the phases of the temporal oscillations of two nearest-neighbours lattice sites 
is a constant a. Thus, there is a phase torsion between the boundaries of the lattice of value r = Net. In 
order to fit the periodic boundary conditions, r = 27tto with m G Z, and, in consequence, a — 2itm/N . 
Hence, we are dealing with an anharmonic plane wave (or a phonobreather) with a wave number equal to a 
that transmits an energy flux along the chain of oscillators [26] . 

Let us remark the similarities of non-time-reversible phonobreathers with q-breathers, which are localized 
solutions in the reciprocal space (see [34] for FPU and [35] for DNLS lattices), q-breathers can be generally 
considered as the superposition of a finite number of normal modes. The solutions considered in the present 
paper, i.e. phonobreathers, are also localized in the reciprocal space, so they can be cast as a special case of 
q-breathers, consisting only of a normal mode characterized by the wavenumbcr a. 

In order to calculate phonobreather solutions we make use of methods based on the anti-continuous 
limit [3, 36], that is, an orbit of frequency cjb for each isolated oscillator is calculated and the coupling 
constant is subsequently varied with a path-following (Newton-Raphson) method. In this paper, we have 
used a Fourier space implementation of the dynamical equations. 

Fourier space methods are based on the fact that the solutions are Tb-periodic (For a detailed explanation 
of these methods, the reader is referred to Refs. [37, 38, 39]). Thus, they can be expressed in terms of a 
truncated Fourier series expansion: 

k m 

u n (t)= ^2 z k ^exp(ik(uj h t + a)) (4) 

k=-k m 

with k m being the maximum of the absolute value of the running index k. In the numerics, k m has been 
chosen as 13. After the introduction of (4), the dynamical equations (1) transform into a set of N x (2k m + 1) 
algebraic equations where the variables are Z = {zk. n }'. 

Fk, n = -ulk 2 z k m + T k [V'(u n )} + e(2z k m - z k ^ x - z k , n +i) = (5) 

Here, T k denotes the Discrete Fourier Transform: 

k 

J r k[V'(u n )]= * +i V'iuitjVexpi-ikubtj), (6) 

j=-k m 
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Figure 1. Position u n (t) and linear momentum p n (t) = w»(t) profiles (left panel) and time evolution (right 
panel) for a phonobreather with k>t = 0.8, £ = 0.05, a = 2tt/3 and a Morse potential. Different curves in 
right panels stand for different lattice sites. Notice that there are only 3 curves in those panels because this 
is the spatial periodicity of the phonobreather. 



where tj is a sample of times that must be chosen equally spaced: 



_ 2?rj 
j ~ wb(2fc m + l)' 



j — km, i • ■ • ) H~^m; 



(7) 



and u n (tj) is calculated from the Fourier coefficients Zf. „ by means of the Inverse Discrete Fourier Transform: 



l n(t 3 ) = ^ Zk, n eyip{ikuj h tj). 



(8) 



k=-k„ 



Fourier space methods provide with an analytical form of the Jacobian J = dF/dZ, whose elements 
are {dFk t n/dzk', n '}- Figs. 1 and 2 show the profiles and time evolution of two examples of phonobreathers 
| with different phases and potentials. 

Phonobreathers can transport energy by means of the phase torsion mechanism, and consequently, there 
is a stationary energy flux as long as a ^ 0, n [7]. With the aid of the Fourier series expansion, the energy 
flux can be expressed as: 



Th Jo 



u n (t)u m (t)dt = -iw h e kzk, n z-k,; 



(9) 



Phonobreathers are characterized by a constant energy density at each lattice site. Thus, it must be 



fulfilled that = Zk. n exp(ifca) with zu,- 

and fulfills: 



and the flux between neighbouring sites is homogeneous 



Jo = Jn.n+i = 2wb£^ k\z k , n \ 2 sin(fca) 



(10) 



fe>i 



f Notice that hereby we will omit, for simplicity, the non-time reversible character of the solutions 
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Figure 2. Position u n (t) and linear momentum p n {i) = ii n {t) profiles (left panels) and time evolution 
(right panel) for a phonobreather with u)^ = 2, e = 0.05, a = 7r/4 and a </> 4 potential. Different curves in 
right panels stand for different lattice sites. Notice that there are only 4 curves in those panels because this 
is the spatial periodicity of the phonobreather. 



2.2. Linear stability equations 

In order to study the linear stability of phonobreathers, we introduce a small perturbation to a given 
solution u n0 of Eq. (1) according to u n = u„q + £ n . Then, the equations satisfied to first order on £ n is: 

Cn + V"( Un0 )£„ + E(2£„ - " Sn-l) = , (11) 

or, in a more compact form: 

AA({u(;)})£ = 0, (12) 

where Af({u(t)}) is known as the Newton operator. In order to study the orbital stability analysis of the 
relevant solution, a Floquet analysis can be performed if there exists 7& € K so that the map u n (0) — > u n (Tb) 
has a fixed point [7]. Then, the stability properties are given by the spectrum of the Floquet operator A^o 
(whose matrix representation is the monodromy) defined as: 

{ ^ (Tb)} ^ =M ( K» ) (13) 



The 2N x 2N monodromy eigenvalues A = exp(k9) are dubbed the Floquet multipliers. This operator 
is symplectic and real, which implies that there is always a pair of multipliers (corresponding to the phase 
and growth modes) at 1 and that the eigenvalues come in duplets {A, 1/A} if they are real or quadruplets 
{A, 1/A, A*, 1/A*} if they are complex. Consequently, if the phonobreather is stable, all the eigenvalues lie 
on the unit circle. 

Equation (12) can be seen as the eigenequation of the Newton operator for the eigenvalue E — 0. Then, 
the eigenequations for the Newton operator are: 

JV({u(t)})C = -> in + V"(u n0 )U + e(2£„ - Cn+i - z„_x) = E£ n (14) 

The Newton operator is periodic in time, and, consequently, its eigenvectors fulfill the Floquet-Bloch 
theorem. This theorem implies that the -E-eigenvalues spread bringing about a set of bands £ v = £ v (0), 
where 9 can be chosen in the First Brioullin zone, i.e. 9 G [— ir, tt]. The set of eigenvalues E v (9) with 9 real 
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is denoted as the v-ih band. The bands are associated to stable solutions as long as 6 V G R. The values of 
9(E) can be obtained by diagonalizing the matrix M. e, which is obtained in a similar fashion to Eq. (13) but 
integrating Eq. (14) for each value of E. The monodromy corresponds obviously to E — and, consequently, 
the Floquet arguments correspond to 9 U (0). Thus, a solution is stable if there are be 2N bands that either 
cross the E — axis or are tangent to it. More details on this theory, called Aubry's band theory can be 
found in Ref. [7] . 



3. Analytical results 

In this section, we will show some analytical predictions about the stability of phonobreathers at low coupling. 
To this end, we start by recalling previous results established by some of the authors for multibreathcrs with 
an arbitrary number of excited sites and continue by applying these results to the prediction of the stability 
properties of phonobreathers with Morse and (j) 4 potentials at low coupling. 



3.1. Previous results 

For the sake of completeness, we recall in this subsection some previous results on the Multibreathcrs Stability 
Theorem (MST) proposed in Ref. [9]. For more details, the reader is also referred to Rcfs. [10, 15]. 

The MST refers to Klein-Gordon systems of the form (1) and estimates the displacement experienced 
by Aubry's bands when the coupling parameter e is switched on. 

Suppose that is a Tb-periodic solution at the anti-continuous limit (e = 0), with p excited oscillators 
and N — p ones at rest (zt° =0). At this limit, there are p degenerated bands tangent to the E = axis at 
(E, 8) — (0, 0). Their curvature is positive (negative) for soft (hard) on-site potentials. The MST can predict 
the displacement of the minimum of this bands {AE}^~q = {E}f~y = s{^}fZo^ with {A} being the set of 
eigenvalues for the perturbation matrix Q. The (non-zero) non-diagonal elements of the p x p perturbation 
matrix Q in reduced form (see below) are defined as 

I fT b /2 

Q n>n±1 = / u° n u° n±1 dt (15) 

M" Anil J-Tt,/2 

with \i n — J T ^ 2 (u^) 2 dt. Only the indexes corresponding to the excited oscillators arc considered and 
they are renumbered them from 1 to p. If we considered all the oscillators, the matrix Q would be in full 
form. Each oscillator at rest adds a row and a column of zeros and, therefore, a zero eigenvalue, which is 
not relevant for the stability properties. 
The diagonal elements are given by 

n,n = • (16) 

With these definitions, we reproduce here the multibreather stability theorem: 



Generalized MST Given a Klein-Gordon system, Eq. (1), a specific multibreather solution at zero coupling 
{it°}, u(t) the corresponding solution at low and positive coupling, {A}^ 1 the eigenvalues of the reduced, 
perturbation matrix Q, with only one zero, then: 
The solution u(t) is stable if: 

a) The on-site potentials are soft and there is not any positive value in {Aj}?~ . 

b) The on-site potentials are hard and there is not any negative value in {\}^Zo ■ 



Stability of non-time-reversible phonobreathers 



7 



We can summarize the stability properties in the following way. If S = 1 means stability and S = — 1 
instability, H = 1 corresponds to a hard on-site potential, and H = — 1 to a soft one, and we define 
sign(Q) = 1 if all the eigenvalues of Q but a zero one are positive and sign(Q) = — 1 if they are negative 
except for the zero one, then 

S = H x sign(Q). (17) 

It there are eigenvalues of different signs, the multibrcather is always unstable. It is important to take into 
account that there is always a zero eigenvalue due to a global phase mode. If there is more than one zero 
eigenvalue, the stability theorem can only predict the instability in the case that there exists at least one 
eigenvalue A that leads to S = —1 in the previous equation (changing sign(Q) for sign(A)), but not the 
stability, as the O-eigenvalue is degenerate. 



3.2. Application to phonobreathers 

We apply hereby the theory recalled in the previous subsection for the system given by Eq. (1). Consequently, 
we must construct the reduced perturbation matrix Q whose elements are given by Eqs. (15) and (16). The 
functions w°(i) that appear in those equations are the solutions of the isolated oscillators submitted to the 
potentials V(u n ), i.e. the solutions of the equations: 

u° n + V'(u° n )=0, (18) 

Let I nm be defined as: 

In,m= I ^ U° n {t)u° m (t)dt, (19) 
J-T b /2 

and the parameters z9„, m as: 

K, m = -7===- (20) 

Substituting Eq. (4) into (19)-(20) and taking into account that the phase difference between 
neighbouring sites a is constant, we get that 



~^^k 2 zl cos(Aia) 



k 

With the aid of the d parameter, the Q matrix can be written as 

Q = VQ {R \ (22) 

where is the Q matrix for a time-reversible phonobreather (i.e. with a = 0) in an homogeneous 
lattice: 

Q ( n R l = 26 n , m -(6 m , n+1 +5 n , n _ 1 ), n = l,...N (23) 

Thus, the eigenvalues of our problem are given by A = i9A# with Xr the eigenvalues of the time-reversible 
lattice. These eigenvalues are given by Xr^ = 4:Sm 2 (mr/N) with n = 0, . . . , N — 1 [9]; consequently: 

A„=4i?sin 2 ^ n = 0,...,N -1. (24) 

Thus, sign(A„) = sign(^). Following Eq. (17), the phonobreather is stable for a soft potential if $ < 
and, for a hard potential, if $ > 0. In the case of a symmetric on-site potential, z 2iU — 0,W <E Z, and 
consequently, ■&(tt/2) = 0. 

We calculate below the expression of $ for a Morse and a hard ^ 4 potential: 
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Figure 3. Stability and instability regions predicted by the MST for Morse (left) and <f> 4 (right) potentials. 
Recall that, for the Morse (0 4 ) potential, instability implies i? < (i9 > 0). 



3.2.1. Morse potential The orbits of Eq. (18) are given by [9] 



u(t) = lo^ 



1 — i/l — CJj^ COS UJb t 



The Fourier coefficients are: 

1 + UJ h 



Zn,n = log • 



(-1)*+* fl-uj h \ k/2 



Thus, we can write: 



with 



J2kLi rk cos (ka) 

1 - CJ b 

r = . 

1 + UJ h 

Using [40], the sums can be performed, leading to: 

^ (1 — r)(cos a — r) 
1 — 2r cos a + r 2 

Thus, 

sign($) = sign(cos a — r), 
In consequence, i9 < (i.e. the phonobreather is stable) if 

1 - w b N 



cos a < r 



1 + UJ h 



(25) 

(26) 

(27) 
(28) 

(29) 
(30) 

(31) 



The main consequence of this condition is that, for a given value of a, there exists a critical value of uj^ 
below which the phonobreather is unstable. This critical value is: 



1 — cos a 
1 + cos a 



= tan 



(32) 
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Consequently, as w b < 1 for oscillators with Morse potential, if a > n/2 phonobreathers are always 
stable, independently of the frequency. 

Fig. 3(left) shows the predictions of the MST for the Morse potential as a function of w b and a. For 
Wb = Wb,c(a), $ = and the degeneracy of A„ cannot be removed; consequently, the MST would not be able 
to predict the (in)stability in that case. 

3.2.2. <p 4 potential The orbit of an oscillator submitted to a 4 hard potential is given by [41]: 



u n {t) = X j T -^cn [ 7f =^, *) = V 7^ cn K J < ( 33 ) 

where cn is a Jacobi elliptic function of modulus n n and K(k) is the complete elliptic integral of the first 
kind defined as: 



r ,z da; 

K(k) = / - AX . (34) 
"'o v 1 — « sin x 

The breather frequency cjb is related to the modulus n through: 

w b = ; • (35) 

2V1 - 2k 2 K(*) V ; 

In order to calculate it is better to use Eq. (20). Taking into account that # depends only on the phase 
difference, we can write: 



7r(l-2/t 2 ) j_ 2K(k) 



dx sn(x)sn(x + q)dn(x)dn(x + a), (36) 

-2K(k) 



with a = 2K(K)a/n. This integral can be evaluated applying [42, identity 171]: 



K 2 sn(x)dn(x)sn(x + q)dn(x + a) = 

— cs(a)ns(a)(l + dn 2 (q)) + cs(q)ns(q)(dn 2 (x) + dn 2 (x + a)) — 

ds(a)(cs 2 (a) + ns 2 (q))[Z(x + a) - Z(x) - Z(a)], (37) 

where Z(x) is the Jacobian elliptic Z function [43]. Then, the value of $ is given by: 

= cs(q)ns(q)[2£(*) - K{n){l + dn 2 (q))] - ds(q)(cs 2 (q) + ns 2 (q))Z(q) 

{-1 + 2k 2 )E(k) + {1-k 2 )K(k) ( ' 

where E(k) is the complete elliptic function of the second kind: 

r /2 i 

E(k) = VI- k 2 sin 2 x dx . (39) 

Jo 

Three relevant values of d correspond to a — tt, a — n/2, a — which lead to $ = 1, $ = and i9 = — 1, 
respectively. This result can be straightforwardly determined from Eq. (21). An analytic determination of 
the sign of $ for other values of a is not possible. However, it is possible to visualize numerically that the 
only zero of #(a) takes place at a = tt/2 being $ < (d > 0) for a < tt/2 (a > n/2). In consequence, 
phonobreathers with a </> 4 hard potential are stable as long as a < tt/2. 

Fig. 3(right) summarizes the predictions of the MST for the <j) 4 potential as a function of uj\, and a. For 
a = tt/2, # = and the degeneracy of A„ cannot be removed, and consequently, the MST could not predict 
(in)stability. 
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4.5- 




Figure 4. Flux between nearest-neighbours for phonobreathers with a = 2ir/3 and a Morse potential. 
Dashed lines represent unstable solutions, for which the flux diverges above (below) a critical value of e 
(«J b )- 



4. Numerical results 



The analytical results of the previous section concern to an infinitesimally small coupling. In this section, 
we explore the existence and stability properties for an arbitrary coupling constant and determine the range 
of validity of the predictions of the multibreather stability theorem. 



4-1. Morse potential 

Before undertaking the existence and stability analysis of the phonobreather, we show the properties of the 
energy flux. Fig. 4 depicts the dependence of Jo with respect to e and LUh for a — 2ir/3; the qualitative 
behaviour displayed in the figure is generic for any value of a. It is observed that for the Morse potential, the 
flux increases (decreases) with the coupling constant (frequency). There is a critical value of e (wt>) above 
(below) which phonobreathers do not exist. For this critical value, the flux diverges. This fact is used in 
Appendix A to calculate the dependence of the critical value of e with Wb ■ 

Additionally, a full stability analysis has been performed by varying e and Wb for a Morse potential and 
choosing two different values of the phase difference a = 2tt/3 and a — Att/9. In the first case, phonobreathers 
are stable for small e and become unstable above a critical value of e. In the second case, phonobreathers 
are unstable for small e as long as Wb > ^b.c! the MST predicts that Wb.c = 0.7041, which fits quite well with 
the numerics; in addition, for large enough e, the phonobreather becomes stable, losing the stability again 
when e is increased further. As explained above, phonobreathers cease to exist above a critical value of s 
because the flux diverges. This critical value of e has been calculated in Eq. (A. 6): 

(40) 



c 4 sin 2 (a/2) 

Fig. 5 depicts the existence and stability ranges for phonobreathers with the above mentioned parameters 
whereas Fig. 6 shows the dependence with respect to e of the arguments and moduli of Floquet multipliers 
for several examples. From these figures, it can be deduced that instabilities and stabilities arise by means of 
cascades of Neimark-Sacker bifurcations (also known as Krein crunches). An analysis of the dynamics shows 
that those instabilities lead to the disappearance of the phase torsion and the time and spatial periodicities 
of the phonobreathers. 
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Figure 5. (s-ojt,) plane showing the existence and stability of phonobreathers with a = 4ir/9 (left) and 
a = 2ir/3 (right) in a Morse potential. Dashed line in left panel corresponds to the MST prediction (i.e. 
stable at the right of the line and unstable at the left) for low coupling. 



We have checked the analytical predictions of MST regarding to Aubry's bands displacement. To this 
end, we show in Fig. 7 the bands for a phonobreather with small e whereas Fig. 8 shows the comparison 
of the predicted vertical displacement given by the MST with the numerical value. We can observe that, 
contrary to time-reversible multibreathers (where the displacement is purely vertical), there is a diagonal 
movement of the bands. In addition, there is a horizontal splitting of the bands that preserves the mirror 
symmetry of the bands with respect to i?-axis. It is also found an excellent agreement with the analytical 
predictions even for high coupling when a — 2ir/3 (the displacement is almost lineal), whereas for a = 47r/9 
there is a serious discrepancy as the bands displacement is non-monotonic and non-linear. 

4-2. Hard 4 potential 

The analysis of the previous case can also be done for the hard 4> 4 potential. We start by showing in Fig. 9 
the properties of the flux for a — tt/A and a — 37r/2; the qualitative behaviour displayed in the figure is 
generic for any value of a. Contrary to the Morse case, the flux has a non-monotonical behavior, increasing 
with e for small coupling and decreasing for large coupling; additionally, the flux is zero for zero coupling 
and for a critical value of e. However, the flux increases with the frequency. As in the Morse case, this 
behaviour allowed us to calculate the critical coupling as a function of the frequency, as shown in Appendix 
A. 

Let us show the results regarding to the stability. Analogously to the time-reversible breathers, there 
exist no stability exchange bifurcations and the stability for any arbitrary e is the same as for small coupling 
except for a — tt/2. Thus, for a > n/2 (a < 7r/2), phonobreathers stability is the same as in the a = n 
(a = 0) case. The only unknown appears at a = tt/2, where the multibreathers stability theorem fails. In 
that case, numerical calculations demonstrate that the stability range depends both on e and cjb- Fig. 10 
displays the dependence of the Floquet multipliers with respect to e for a — 2ir/3, a = tt/2 and a — tt/A. 
In all cases, for a given u>b, there is a critical value of e where the continuation fails. At this value the 
flux vanishes, contrary to the Morse potential, for which the flux diverges. This critical value of e has been 
calculated in Eq. A. 5: 




Figure 6. Dependence of the Floquet multipliers with respect to e for phonobrcathers in a Morse potential. 
Parameters are: m = 0.68 and a = 47r/9 (top panels), Ub = 0.8 and a = An/9 (middle panels) and Ub = 0.8 
and a = 2tt/3 (bottom panels). 
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Figure 8. Dependence of the location of the minima of the two bands closest to E = (represented by 
E\ and E2) with respect to e for phonobreathers in a Morse potential with a = 4n/9 (left) and a = 2ir/3 
(right). In both cases, u) h = 0.8. Dashed lines correspond to the MST predictions and full lines to the 
numerical values. Solutions at left, where the discrepancy is high, are unstable for small but finite coupling. 



Fig. 11 shows the existence range for those values of a. It can be deduced that the existence range is 
enlarged when a decreases. Besides, for a = tt/2 it is observed that the stability properties depends on the 
frequency and coupling and that phonobreathers are always unstable for cjb < 3. 

Fig. 12 and Fig. 13 show, respectively, examples of bands and the MST predictions, for a — 2tt/2> 
and a = 7r/4. There is an excellent agreement of the numerical results with the analytical predictions. 
Additionally, Fig. 14 focuses on the a = tt/2 case with Wb = 2. From the last figure, we can observe that the 
first and second minima of the bands adjust respectively to the parabolas — 0. 015509s 2 and — 0.056476e 2 . This 
fact confirms that a second order perturbation theory is needed in order to predict the stability properties 
for a = tt/2. 




Figure 9. Flux between nearest-neighbours for phonobrcathers in a rf> 4 potential with a = n/A (top panels) 
and a = 3ir/2 (bottom panels). Dashed lines represent to unstable solutions. Note that the flux tends to 
zero when the continuation in e fails. 



5. Conclusions 

The study carried out in this paper has been twofold: on the one hand, analytical predictions on the 
stability of non-time reversible phonobreathers has been established by making use of the multibreathers 
stability theorem; on the other hand, a numerical study on the existence and stability properties of those 
phonobrcathers has been performed. 

The multibreathers stability theorem predicts that there is a critical value of the phase difference in 
phonobreathers that separates stable and unstable states. This critical value depends on the phonobreather 
frequency for lattices with harmonic coupling and Morse on-site potential whereas, for the <p A potential, the 
critical value is 7r/2 and does not depend on the frequency. At this particular value of the phase, the theorem 
cannot make any prediction about the (in) stability. 

The analytical findings have been confirmed by means of numerical analysis, showing in addition that 
the stability changes above a given value of the coupling constant. It has also been confirmed that, in most 
cases, the Aubry's band displacement predicted by the multibreather stability theorem for small coupling, 
is valid even for relatively high coupling. Apart from this, it has been observed that phonobreathers cease 




Figure 10. Dependence of the Floquet multipliers with respect to e for phonobreathers in a ifi 4 potential. 
Parameters are: aj^ = 2 and a = 2tt/3 (top panels), u)-^ = 2 and a = n/2 (middle panels) and uj^ = 2 and 
a = n/A (bottom panels). 





Figure 12. Bands for a phonobreather with e = 0.05, wt, = 2 and a = 2ir/3 (left) and a = n/4 (right) and 
a </> 4 potential. The number of particles is N = 21 (left) and N = 24 (right). 



to exist above a critical value of the coupling. This value coincides with a flux divergence for the Morse 
potential or a flux vanishing for the hard </> 4 potential. Finally, it has been found that for the </> 4 potential 
with a = 7r/2, the (in)stability depends on the frequency even for infinitcsimally small coupling, being 
additionally confirmed that a second-order perturbation theory is needed in order to predict the (in) stability 
at this particular value of the phase. 

The above mentioned results open new perspectives for future work. For instance, we have started to 
study the properties of vortices in square lattices, which can be pictured as multibreathers with a = tt/2. It 
can be interesting to apply the multibreathers stability theorem to two structures considered by Cretegny and 
Aubry: (i) breather "rivers" (percolating cluster) firstly and (ii) phonobreather with vacancy sites. Finally, 
we have also started to study the effect of long-range interactions in the stability of non-time reversible 
structures. 
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Figure 13. Dependence of the location of the maxima of the two bands closest to E = (represented by 
Ei and E2) with respect to e for phonobreathers in a </> 4 with a = 2-7r/3, N = 21 (left) and a = 7r/4, TV = 24 
(right). In both cases, cjb = 0.8. Dashed lines correspond to the MST predictions and full lines to the 
numerical values. Dashed lines correspond to the MST predictions and full lines to the numerical values. 




Figure 14. Phonobreather in a </> 4 potential with o^, = 2, a = n/2 and TV = 20. (Left) Band diagram for 
£ = 0.05. (Right) Dependence of the location of the maxima of the two bands closest to E = with respect 
to e. 
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Appendix A. Calculation of the critical coupling 

In this appendix we calculate the critical value of the coupling parameter e for the on-site potentials 
considered in the paper, i.e. the Morse and the (j) 4 potentials, by making use of the rotating wave 
approximation (RWA) (see e.g. [44]). This method consists in considering all the Fourier coefficients zero 
but the first one, as this is usually considerably larger that the other coefficients. The RWA gives often 
surprisingly good results, taken in account its simplicity. 
Let us consider equation Eq. (1): 

u n + V'(u n ) + e(2u n - u n+ i - u n -i) = 0, n = l,...,N. (A.l) 

We can write u n and V'{u n ) in terms of their fourier coefficients z k _ n and k as 

k m k m 

u n (t)= z k, n e*p(ikuj h t); V'(u n (t))= ^ V k\n cx P( ikuj bt), (A.2) 

As Zk, n ±i — Zk, n exp(±ia), substituting in Eq. (A.l) we obtain 

-k 2 uol z k , n + V^ n + Ae z Kn sin 2 (a/2) . (A.3) 

Let us calculate the coefficients V kn for the generic 4 potential V(u n ) = \u 2 n + jsu^. Then 
V'(u n ) =u n + sul or 

V'(u n ) = ^ z k , n exp(ikuj h t) + s( ^ z fe ,„ exp(ifcw b i)) . (A.4) 

k——k m k——k m 

As u n (t) is real zi jTl = z*_ l n and, therefore, within the RWA u n (t) = z\ yTl exp(iubt) + %\ n exp(— iuj^t) = 
2\zi tn \ cos(uJbt + 4>i, n ) with 4>\ t n = arg(zi.„). For a given n = p, we can always set 4>i,p = by choosing the 
appropriate time origin; that is z\ tP e M and z\ tP = 2_i jP . We only need to find the coefficients for the site 
n = p as we can deduce any other one by using z k , n +p = z k ^ p exp(ma). Let us write u(t) = u p (t), z k = z ktP 
and so on, for simplicity. Then 



V ' (u(t)) = [z\ exp(iwb^) + z\ exp(— iujht)) + s(z\ exp(iwb^) + z\ exp(— iuj^t)) 3 = 
Zi(exp(iuibt) + exp(-iu)bt)) + szf(exp(iu>bt) + exp(— ioj^t)) 3 = 

Zi(exp(iujbt) + oxp(— iuobt)) + szl(3[exp(iujbt) + exp(— iu>bt)] + [exp(3iu;bi) + exp(— 3iu)bt)}) = 
(zi + 3s z 2 )(exp(iu>bt) + exp(-iiObt)) + szf(exp(3iu}bt) + cxp(— 3iujbt)) . 

Therefore, the coefficient of exp(iu)bt) is V{ = (z\ + 3s z\). Substitution in Eq. (A.3) for n = p 
(omitted) gives — uj 2 zi + (z\ + 3s z\) + Aezi sin 2 (a/2) = 0. This equation has two solutions z\ = and 
z\ = (—i^b + 1 + 4esin 2 (a/2)), which becomes also zero for the critical value of the coupling parameter e 
given by: 

LO 2 - 1 

£ c = : .\, m ■ (^potential) (A.5) 
4 sin (a/2) 

Therefore, z\ tP = z\ = and so are z\, n and u n within the RWA. The flux Jo also becomes zero, according 
to Eq. (10). This result does not depends on the value of s, which measures dc degree of nonlinearity of the 
potential. As any hard symmetric potential can be approximated by the cf> 4 potential with some s > 0, e c is 
a good approximation for all of them. However, for a soft symmetric potential s is negative, u>b < L an d as 
£ > there is no value of the coupling that nullifies z\ and the flux cannot be nullified. 
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For the Morse potential, the RWA approximations is not so useful. From the numerical results we 
know that the critical coupling s corresponds to a divergence of the flux J , and therefore to a divergence 
of u n and its first fourier coefficient Then, V'(u n ) = 2(exp(— 2u n ) — exp(— u n )) tends to zero when 

u n — > oo and so does V{ n . By substitution in Eq. (A. 3) for k = 1, neglecting V{ n and simplifying zi iTl we 
get — lu^ + 4e sin 2 (a/2) = and the critical coupling is 

, ,2 

(Morse potential) (A. 6) 



C 4 sin 2 (a/2) ' 

We have not been able to deduce the reason why u n becomes infinity. 
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